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A directed polymer is considered on a flat substrate with 
randomly located parallel ridges. It prefers to lie inside wide 
regions between the ridges. When the transverse! width 
W = exp(AZ/^''^) is exponential in the longitudinal length 
L, there can be a large number ~ expL^^^ of available wide 
states. This "complexity" causes a phase transition from a 
high temperature phase where the polymer lies in the widest 
lane, to a glassy low temperature phase where it lies in one 
of many narrower lanes. 

Starting from a uniform initial distribution of independent 
polymers, equilibration up to some exponential time scale in- 
duces a sharp dynamical transition. When the temperature is 
slowly increased with time, this occurs at a tunable tempera- 
ture. There is an asymmetry between cooling and heating. 

The structure of phase space in the low temperature non- 
equilibrium glassy phase is of a one-level tree. 
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The glass transition is caused by the appearance of a 
multitude of long lived states, which prevents exploration 
of the whole phase space. These effects are so strong 
that in practice one can only observe precursor effects. 
Experimentally one observes a dynamical freezing around 
a tunable temperature. 

The ergodic theorem says that time-averages may be 
replaced by ensemble averages. It is widely believed that 
the inherent dynamical nature of the glass transition im- 
plies that there is neither need nor chance for a thermo- 
dynamic explanation. However, since so many decades in 
time are involved, this is an unsatisfactory point of view. 
Jackie and Palmer |^ have therefore considered a glass 
as a system with many states ("components") a, each 
having its own component free energy Fa(T). Within a 
state thermodynamic equilibrium is established on mod- 
erate time scales, while the system may escape to other 
states on much longer times. It is often assumed that the 
states occur with Gibbs weight Pa = exp{— PFa{T)) / Z . 
From these assumptions it was recently pointed out that 
the sharp dynamical glassy transition of a typical mean 
field model, the p-spin interaction spherical spin glass, 
can be explained thermodynamically: this transition 
is driven by the complexity or configurational entropy 
(i.e. the logarithm of the number of relevant states). 

Here we wish to study this picture of the glassy transi- 
tion in a simple, exactly solvable model without frustra- 
tion, both at the static and the dynamical level. 

Consider a directed polymer (or an interface without 



overhangs) z{x) which lies in the region 1 < x < L and 
1 < z < W of the square lattice with unit lattice param- 
eter. In the Restricted Solid-on-Solid approximation the 
interface can locally be flat {z{x + 1) = z{x); no energy 
cost) or make a single step {z{x + 1) — z{x) — ±1; en- 
ergy cost J). The partition sum of this system, subject 
to periodic boundary conditions, can be expressed in the 
W X W matrix T that transfers the system from x to 
x + 1 

w 

Z = Tre-^^ ^TrT^ = ^ (1) 

w — l 

where A^u is one of the eigenvalues of T. For a pure sys- 
tem at temperature T ~ 1/(3 Fourier analysis gives the 
eigenvalues A{k) = 1-1- 2e~'''^cosfc. The largest eigen- 
values occur at small momentum, A(fc) « exp[— — 
rfc^/(27r^)], which deflnes the bulk free energy density 
fsiT) = -Tln(l + 2e-'^-') and the stiffness coefficient 
r(r) = 2n^e-f^J /{l + 2e-~^J). 

As depicted in Figure 1, we shall consider the situ- 
ation of randomly located potential barriers parallel to 
the a;-axis, so V{x,z) = V{z). Hereto we assume bi- 
nary disorder, V{z) — with probability p = exp(— /i) 
or V{z) = Vi > with probabihty 1 - p. Eq. (|l|) is 
dominated by the largest eigenvalues. It is well known 
that they occur due to Lifshitz-Griffiths singularities. 
[|| 11 §1 They are lanes of width £ > 1 in which all 
V{z) = 0, bordered by regions with V{z) ^ 0. These 
dominant configurations are the "states" or "compo- 
nents" of our system. (In spin glass theory such states 
are called "TAP-states" .) The eigenf unction centered 
around Za has the approximate form cos[7r(2; — Za)/ia] 
inside the lane, while outside the lane it decays essen- 
tially exponentially due to the disorder. These states 
can thus be labeled by o = {za,ia)- Since k n/£a their 
free energy follows as Fa = Fg^ where 

PFe^-LlnAi-)^PfBL+— (2) 

The number of regions with £ successive sites with V = 
can be estimated by their ensemble average, A/i — W{1 — 
pYp^. Since the relevant £'s will turn out to be of or- 
der L^/'^, the interesting situation occurs when we choose 
W = exp(AL^/^). The states with width i have a config- 
urational entropy or complexity Xg = In A/f ~ \L^/^ — fii, 
where fi = — In p. For large L we may restrict the par- 
tition sum to these dominant states. We thus evaluate, 
instead of eq. (j^), the "TAP" partition sum 



1 



Z : 



The total free energy 



(3) 



= - In Z = L/3/b + 2^ - AL^/^ + iii (4) 

has to be optimized in I. The largest i which occurs in 
the system can be estimated by setting Mt ~ 1, yielding 



(5) 



It is a geometrical length, independent of T. Let us in- 
troduce 7(r) = (r(r)/i2)i/3/A. The free energy of this 
widest state then reads 



f3F = LPfB + ^XL'/'-f{Tf 



(6) 



It follows from M) that at low enough T the optimal 



idth is smaller than 



max ; 

rixi/3 



7(^)4 



The free energy of this phase is 

/3i^ = L/3/B + iALi/3(37(T)-2) 



(7) 



(8) 



In the temperature interval where 7(T) > 1 the interface 
lies in the non-degenerate widest lane. For 7(T) < 1 it 
lies in one of the N't* ^ 1 optimal states, each of which 
has a higher free energy than the widest lane; this free en- 
ergy loss is more than compensated by their complexity. 
The system thus undergoes a static glassy transition at 
temperature Tr- — J/ ln(27r^^^A^^-2), where 7 = 1. For 
any finite L there is also a very low temperature regime 
T < 1/logL, where the interface is essentialy straight 
and can lie anywhere in the system, and a high temper- 
ature regime T > logL, where the potential barriers are 
ineffective, and the interface shape is truly random. 

When considered as function of T = 1/T this situ- 
ation is reminiscent of the p-spin interaction spin-glass 
and to the Random Energy Model |Q. The very- high 
T regime extends up to T ~ In L, in which regime a 
gradual freezing takes place in TAP states of width £* 
much larger than unity. This smeared transition is re- 
lated to the sharp dynamical transition at some tem- 
perature Ta > Tk of mean field spin-glass models. In 
the regime Tk < T ^ InL the system is frozen in 
TAP states of appropriate degeneracy; a similar static 
phenomenon was found for spin-glass models in the 
regime Tk < T < Ta- Like in these models, below the 
"Kauzmann" temperature Tk there is only an essentially 
non-degenerate state. This a manifestation of the "en- 
tropy crisis" of glasses and glassy systems. 



The internal energy of a state of width £ is 

At the transition point 7 = 1 it is continuous, sim- 
ply because £* £max- The free energy is a ther- 
modynamic potential, and yields the same value for U . 
There also holds the relation dJ7/dT = TdS'/dT, where 
S = St,(T) +I{T) = -dF/dT is the total entropy On 
the side 7 < 1 the free energy (^) deviates quadratically 
from (ra), leading to a higher specific heat. It reads 



^ = = + ^^^^^^ + ^Ail/3(j.5^^)2 (10) 

It exceeds the component averaged specific heat of this 
phase, C = Y^aPaCt^ = LcB + {L/2e)dTdpT. The spe- 
cific heat is larger in the glassy phase than in the param- 
agnet, because the free energy is lower. When considered 
as a function T, the specific heat makes a downward jump 
on cooling through Tr-, as occurs in realistic glasses. 

When monitoring the internal energy of one polymer 
as a function of time, as is easily done in a numerical 
experiment, one obtains essentially a noisy telegraph sig- 
nal. Each plateau describes trapping of the polymer in 
one lane for some definite time. The variance of the noise 
in the internal energy on this plateau is equal to T^C, 
implying fluctuations of order L^^^. From time to time 
the polymer moves to another lane, causing additional 
noise. The variance of the total noise equals T^C, and it 
indeed exceeds C by an amount of order L^/^ . 

We now consider dynamics. On appropriate time 
scales, our system can be viewed as a set of deep states 
(traps) a located in lanes of width £a around centers z^, 
with associated free energies Fa = Fg^ given in eq. (2). 
These minima are separated by very wide regions (sep- 
aration ~ expL^/"^) with fully random potential, that 
builds a barrier. Between traps a and a + 1 there is a 
free energy barrier determined by the intermediate state 
of highest free energy. Let us call its free energy i?a; it 
will typically lie at distance L^^^ below the maximal free 
energy Lf^ax = -LTlnA™„ = -LTln(l - 2e-'''') [|. 
The free energy barrier for the polymer to move from 
state a to state a -|- 1 is thus Ba — Fg^ , while the barrier 
for moving to the left is Ba-i — Fg^. For a statistical de- 
scription of dynamics, we assume that many independent 
polymers are present, of which the units make random 
thermally activated moves. On appropriate time scales 
one then gets a master equation for the probability Pa (t) 
that the center of a polymer is inside the a'th state. 



Here to is the attempt time for a move of one polymer 
segment. This model is a combination of the random 
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jump-rate model and the random bond models studied 
before. [ p^|JTl] ], see for a review. The stationary state 
is independent of the barriers Ba '. 



-I3F. 



v7 = 



(12) 



The denominator Zpq is equal to the thermal TAP parti- 
tion sum. So our master equation exactly reproduces the 
Gibbs distribution discussed before. 

Let us now consider the motion of polymers in an equi- 
librium ensemble, after making some non-essential sim- 
plifications of the system. First we slightly modify the 
actual height of the barriers by setting Ba = Lfmax, 
thus neglecting their L^/^ deviations. Next we assume, 
for some fixed Iq (l <^ io ^ L^^^), that all relevant 
deep traps (and further a lot of shallow traps) are lo- 
cated at positions that are multiples of We^, = exp/x^Oi 
and we only consider those states. Their number is 
TVfo = exp(AL^/'^ — /i^o)- Using a result of Haus et al, we 
know that diffusive behavior is exact at all times, pi 



{6zm ^ {{z{t) - z{o)r) = Dt 

involving the diffusion coefhcient 
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(13) 
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with F given in eq. (El). As the barriers have a height that 
deviates by order from a fixed value L{fmax — fs), 
we introduce the logarithmic time variable t [p3 



t{r) 



(15) 



T ranges from ~ —L'^^^, where t ^ to, up to 0(1), where 
t ~ io exp(L -I- L^/^) and the interesting physics occurs. 

Let us now consider T < Tk. At given r the 
polymers have a time ^(t) to make moves. The typ- 
ical deviation follows from (|l^) and(p^) as logj^zj ~ 
\L^'^{t + ijimax + 7^^^„^/2F)/2. This we can relate 
to the distance between optimal states Wi = expfii and 
the number of them expli„ter(r), reached at time ^(t): 
\Sz\ ~ WiexpXinter- It defines the intercluster complex- 
ity, i.e. the configurational entropy of the states reached. 



T (t\ — i'- -I- ' max 
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(16) 



The equilibrium relation I = 'y(T)imax tells us that 
for r = 7/2 the next optimal state will typically have 
been reached, so that the population of the states is in 
accord with the equilibrium prediction, and thermody- 
namic equilibrium is achieved. We consider the bunch of 
exp2i„ter(T) states as one cluster. There are eiqylintra 
of these clusters, defining the intra-cluster complexity 

2:^„t™(r) =I-2:™ter(T) = ALi/^(l- ^7- ir) (17) 



As time progresses, more and more optimal states will 
be visited, expressed in the growth of Imter ■ This natu- 
rally occurs by a decrease of the number of independent 
clusters expli„fra- At the ergodic timescale set by 
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T>Tk 



(18) 



ergodic behavior has been reached. Each polymer has 
been able to visit the whole system. The system is 
one big cluster {lintra = 0) that contains all relevant 
states (Xinter — I)- Summarizing, so far we have found 
the regime B (0 < T < T^; l{T)/2 < r < T,rg{T)), 
where more and more optimal states are visited, and the 
regimes C (T < Ta-; r > Te^g(T)) and D ( T > T^; 
T > Terg{T)) of crgodic behavior. Thermodynamic equi- 
librium is already reached in regime B, before the system 
is ergodic. Thermodynamic relations involving deriva- 
tives (e.g. dF/dT, dU/dT, dS/dT) are satisfied. 

We now consider dynamics in the initial time regime 
A (r < 7(r)/2 for T < and r < TergiT) for T > Tk). 
Starting from a uniform initial distribution, where many 
narrow states are populated, the polymers will move to 
wider and wider states. At times t{T) with t < 7/2 a 
polymer will visit a region of width dz{T) where it picks 
up the broadest state available. The non-degeneracy of 
this state (TinteriT) — 0) yields, now using eq. ( p^ ) as an 
estimate, the best width reached up to then. The com- 
plexity of these clusters is Tintra — \L^^^ — Ai^(T). This 
behavior may be expressed in terms of a dynamical par- 
tition sum at time scale ^(t), where the system is split up 
in independent clusters c = 1, ■ ■ ■ , A/'f(r) of width W^^T-y. 

Z{t) = E^'r' Zc{Wtir)) where each of the Zc(VF)'s is 
as in eq. (1). This results in the dynamical free energy 
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Let us consider T < Tk ■ The dynamical free energy has a 
minimum, which is approached in the limit r — > "f(T)/2. 
The value of the minimum coincides with the static value 
(P) . At that timescale a polymer will typically have found 
a state of optimal width £{t) = jimax- The behavior at 
larger times, r > 7/2, was already discussed above, and 
matches nicely with this short time dynamics. 

In order to compare with cooling experiments in realis- 
tic glasses, we must consider a heating experiment (which 
is a cooling experiment in the variable T). The temper- 
ature changes slowly with time, T = T{t). It defines the 
inverse function t{T), that characterizes the heating tra- 
jectory. Due to the L-dependence in (p^), r will start at 
^ — i^2/3 -fQj, g]-Qall t, but when r = 0(1), it need not be 
a monotonically increasing function of T. Approaching 
Tk from below under appropriate conditions, a dynam- 
ical transition will occur at temperature Tf < Tk (it is 
a freezing transition in terms of T). This temperature 
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is set by t(T/) = 7(T/)/2, where, starting from small 
widths, the dynamically achieved width equals the opti- 
mal width, viz. £{T{Tf)) — ^{Tf)£max- ^^Tf the internal 
energy is continuous. For T < Tf the specific heat reads 

C = csL + XL^/^ (^^drd,, + 3:^drln^^ (20) 

For T{t) > Tf the width £ will be equal to the optimal 
value j{T{t))imax and the last term vanishes. In the 
initial time regime, where ^ < 7, this term is positive. 
At Tf the specific heat thus performs a jump, with the 
higher plateau value on the low temperature side, as oc- 
curs for the static behavior at Tk- The height of this 
jump depends on the heating rate, and it vanishes for 
small enough rate. For large L^^^ the jump is sharp in 
the logarithmic time variable r. 

In the dynamical regime A thermodynamic relations 
involving addition, such as F — U — TS remain valid, as 
they still arise from saddle point analysis. But F is not 
a thermodynamic potential. Derivative relations, such as 
d/3F/d/3 = U and dS/dT = C/T are violated. 

There is no symmetry between cooling and heating ex- 
periments. In order to have a decreasing T{t) eq. ( |l5| ) 
tells us that r must be of order —Lp'/'^, whereas it is typ- 
ically of order unity for heating. Equating ( p^ to zero 
then leads to widths £ ^ 0{L'^), thus completely reinitial- 
izing the relaxation. Such a phenomenon was observed 
upon heating in spin glasses and explained in terms of 
hierarchy of phase space. 

Our dynamical analysis puts forward the picture that 
the hierarchical structure of phase space, here having 
the structure of a one level tree and reminiscent of one 
step replica symmetry breaking, is a dynamical effect. 
At given timescale only nearby states can be reached, so 
there are many clusters, each having one state (regime A) 
or many states (regime B). At larger times other states 
can be reached, thus leading to larger clusters. For times 
larger than the ergodic time of a large but finite system 
(regimes C and D), all states are within reach, i.e. there 
is one cluster only, containing many states (regime C) or 
one state (regime D). In the thermodynamic limit L — > 00 
before t — > 00, phase space splits up in truly disjoint sets. 

In contrast to most realistic glasses, our model has 
no crystal state. Further study of the dynamics should 
show whether stretched exponential decay and the Vogel- 
Fulcher law for the relaxation time occur. 
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FIG. 1. A directed polymer can move on a substrate with 
parallel potential barriers. It prefers to lie in wide lanes be- 
tween the barriers. 
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FIG. 2. Dynamical phase diagram for A = = J = 1. 
r = ci + C2 In t is logarithmic in time, T is temperature and 

Tk the static Kauzmann transition temperature. A: initial 
regime; B: glassy regime with many states at thermodynamic 
equilibrium; C: ergodic glassy regime with many states; D: 
ergodic regime with one state. When heating the system from 
short times (large negative r) and small T towards larger t 
and T, a dynamical (or "experimental") transition occurs at 
T — Tg where the A-B or A-D boundary is crossed. Crossing 
the C-D boundary is only possible in the ergodic regime; it 
induces a static "Kauzmann" -transition. 



